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Abstract 

We present a polyhedral algorithm to manipulate positive dimensional solution sets. Using 
facet normals to Newton polytopes as pretropisms, we focus on the first two terms of a Puiseux 
series expansion. The leading powers of the series are computed via the tropical prevariety. 
This polyhedral algorithm is well suited for exploitation of symmetry, when it arises in systems 
of polynomials. Initial form systems with pretropisms in the same group orbit are solved only 
once, allowing for a systematic filtration of redundant data. Computations with cddlib and 
Sage are illustrated on cyclic n-roots polynomial systems. 

1 Introduction 

We consider a polynomial system /(x) = 0, x = (xi, x 2 , . . . , x n ), f = f 2 , . . . , f n ), fi 6 C[x], 
i = 1,2, ... ,n. While in many applications, the coefficients of the polynomials are rational 
numbers, we allow the input system to have approximate complex numbers as coefficients. As 
/(x) = has as many equations as variables, we expect in general to find only isolated solutions. 
In this paper we focus on systems with particular choices of the coefficients so /(x) = has a 
space curve as solution. 

Our approach is based on the following observation: if the solution set of /(x) = has a 
space curve, then this space curve extends from C* = C \ {0} to infinity. In particular: the space 
curve intersects hyperplanes at infinity in isolated points. We start our series development of 
the space curve at these isolated points. We applied this approach to find common factors of 
two polynomials in [1] and a first attempt to generalize this approach to any dimension occurred 
in [23]. Our approach can be interpreted as a symbolic- numeric version of the effective proof of 
the fundamental theorem of tropical algebraic geometry in |12j . 



'This material is based upon work supported by the National Science Foundation under Grant No. 0713018 and 
Grant No. 1115777. 
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The hyperplanes at infinity are associated to the facets of the Newton polytopes of the poly- 
nomial system. In [TO], Newton polytopes are introduced as compactifications of logarithms of 
algebraic sets (also known as amoebas). Formally we denote a polynomial p £ C[x] as 

p(x) = ^c a x a , c a GC\ y* = x?x% (1) 

and we call the set A of exponents the support of p. The convex hull of A is then the Newton 
polytope of p. To every facet of P corresponds a hyperplane at infinity and for every facet we 
can then introduce projective coordinates [22] . 

In this paper we will make various significant assumptions. First we assume that the space 
curves we consider are reduced, that is: free of multiplicities. Moreover, the space curve is in 
general position with respect to the first coordinate plane. More precisely: we assume that the 
space curve is not contained in the coordinate plane defined by X\ = 0. Thirdly, we assume the 
space curve intersects the plane x\ = at regular solutions. The special role of x\ is reflected in 
the normal form of Puiseux series: 

I Xl = ^ (2) 
\ x i = t v >(y i + z l t w >), i = 2,...,n. { ' 

where the leading powers v = (v\, V2, ■ ■ ■ , v n ) equal to what is called a tropism. The leading 
coefficients yi G C* are defined as solutions of the initial form system defined by the tropism v. 

Our approach consists of two stages: 

1. The computation of the pretropisms, which are candidates for the leading powers of the 
Puiseux series [14], [15], [25]. We describe this computation in section 4. 

2. The computation of the leading coefficient and the second term of the Puiseux series, de- 
scribed in the fifth section. 

The second term of the Puiseux series indicates the existence of a space curve. The focus of our 
approach is on systems that have (at least) as many equations as variables, i.e.: for system where 
one would normally not expect to have space curves as solutions. 

An obvious advantage of our approach is the exploitation of symmetry. If the system is 
invariant to permutation of the variables, then it suffices to compute only the generators of the 
solution orbits. We then develop the Puiseux series only at the generators. As we illustrate first 
on the cyclic 5-roots system, this approach of exploiting symmetry is not limited to systems with 
space curves as solutions. 

Our running example is one family of polynomial systems. While still preliminary and limited 
in scope, we are encouraged by the results on the cyclic 12-roots problem. In particular, we obtain 
exact representations for the solution curves derived from numerical homotopy solvers. Our results 
on cyclic 12-roots correspond to [16J. 

2 the Cyclic n-roots Problem 

Our running example throughout this paper is the cyclic n-roots problem, an academic benchmark 
for polynomial system solvers, see e.g: [3], [6], [8], [13] . For n = 3, the system originates naturally 
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from the elementary symmetric functions in the roots of a cubic polynomial. For n = 4, the 
system is more interesting: 



/(*) 



X\ + X 2 + X 3 + ^4 = 
XiX 2 + X 2 X 3 + X3X4 + X4X1 — u 
X1X2X3 + X 2 X 3 X 4 + X3X4X1 + X4X1X2 

X1X2X3X4 — 1 = 



= 



= 



(3) 



We illustrate the application of our approach to this system. Terminology and precise defi- 
nitions are in later sections. There is only one tropism v = (+1, —1, +1, —1) defining the initial 
form in v (/)(z) = 0: 



m. 



•(/)(x) 



x 2 + x 4 = 
xix 2 + x 2 x 3 + X3X4 + X4X1 
X2X3X4 + X4X1X2 = 
X1X2X3X4 — 1 = 



xi 

X3 
X4 



y^y* 



(4) 



The system in v (/)(y) = has two solutions. These two solutions are the leading coefficients in 
the Puiseux series. In this case, the leading term of the series vanishes entirely at the system so 
we write two solution curves as (t, — — t, and (t, — t, — i -1 ) . To compute the degree 
of the two solution curves, we take a random hyperplane in C 4 : cq + c\X\ + C2X2 + C3X3 + C4X4 = 0, 
Ci G C*. Then the number of points on the curve and on the random hyperplane equals the 
degree of the curve. Substituting the representations we obtained for the curves into the random 
hyperplanes gives a quadratic polynomial in t (after clearing the denominator so there are 
two quadric curves of cyclic 4-roots. 



3 Exploiting Symmetry 

We illustrate the exploitation of permutation symmetry on the cyclic 5-roots system. All Newton 
polytopes of the cyclic 5-roots system are in general position. Therefore, the mixed volume is 
sharp and its value is 70. We, hence, know that the solution set for the cyclic 5-roots systems 
consists of 70 isolated solutions. Adjusting polyhedral homotopies to exploit the permutation 
symmetry for this system was presented in |24| . Here we present an alternative approach. 

If we consider the first four equations of the cyclic 5-roots system Cs(x) = 0, then we have 
solution curves. Consider the first four equations of Cs(x) = 0: 

x + xi + x 2 + x 3 + x 4 = 

x xi + x x 4 + xix 2 + x 2 x 3 -I- X3X4 = 
\ ^ (5) 

x xix 2 + x xix 4 + X0X3X4 + X1X2X3 + X2X3X4 = 

k x xix 2 X3 + x xix 2 X4 + X0X1X3X4 + X0X2X3X4 + X1X2X3X4 = 

where v = (l,l,l,l,l). As the first four equations of C5 are homogeneous, the first four equations 
of C5 coincide with the first four equations of in v (C5)(x) = 0. Because these four equations are 
homogeneous, we have lines of solutions. After computing representations for the solution lines, 
we find the solutions to the original cyclic 5-roots problem intersecting the solution lines with the 
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hypersurface defined by the last equation. In this intersection, the exploitation of the symmetry 
is straightforward. 

The unimodular matrix with v = (1, 1, 1, 1, 1) is given by 



M 



11111 
10 
10 
1 
1 



(6) 



The corresponding unimodular coordinate transformation, denoted as x = z , is defined as 

x = ZQ, xi = zqz 1} x 2 = z z 2 , x 3 = z z 3l x 4 = Z Z4 



(7) 



Applying x = z M to the initial form system ([5]) gives 
in v (C 8 )(x = z A/ ) = 



Z\ + Z 2 + Z 3 + Z4 + 1 = 
Z\Z2 + Z 2 Z 3 + Z 3 Z 4 + Z\ + Z4 = 
Z\Z2Z 3 + Z2Z3Z4 + Z\Z2 + Z1Z4 + Z3Z4 = 
ZiZ 2 Z 3 Z 4 + ZlZ 2 Z 3 + Z1Z2Z4 + Z1Z3Z4 + Z 2 Z 3 Z 4 



(8) 



The system has 14 isolated solutions of the form z\ 
we let zq = t, in the original coordinates we have 



Cl, Z 2 = C 2 , Z 3 = C 3 , Z4 = c 4 . 



If 



Xq = t, X\= tCl, X2 = tC2, X 3 = tC 3 , X4 = tC4 



0) 



as representations for the 14 solution lines. 

Substituting Q into the omitted equation xox±X2X 3 X4 — 1 = 0, yields a univariate polynomial 
in t of the form kt 5 — 1 = 0, where k is a constant. 

Among the 14 solutions, 10 are of the form i 5 — 1. They account for 10x5 = 50 solutions. There 
are two solutions of the form (— 122.99186938124345)t 5 — 1, accounting for 2 x 5 = 10 solutions 
and an additional two solutions are of the form (— 0.0081306187557833118)t 5 — 1 accounting for 
2 x 5 = 10 remaining solutions. The total number of solutions is 70, as indicated by the mixed 
volume computation. 

Existence of additional symmetry, which can be exploited, can be seen in the relationship be- 

-0.0081306187557833118. 



tween the coefficients of the quintic polynomial, i.e. (_ 122 99186938124345) 

The next cases in the family that admit positive dimensional solution sets are for n = 8 and 
n = 9 as both 8 and 9 has a square divisor [2]. With more efficient implementations of the 
approach described in this paper, also the cyclic 12-roots problem comes within reach. 

That the first n — 1 equations of cyclic n-roots system give explicit solution lines is exceptional. 
In general, we plan to use the leading term of the Puiseux series to compute witness sets for the 
space curves defined by the first n — 1 equations. Then via the diagonal homotopy [T7] we can 
intersect the space curves with the rest of the system. While the direct exploitation of symmetry 
with witness sets is not possible, with the Puiseux series we can pick out the generating space 
curves. 
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4 Computing Pretropisms 



Let v/0 and denote (a, v) = a\v\ + a 2 v 2 + • • • + a n v n . Then in v (/) is the initial form of / in 
the direction of v: 

in v (/)= Yl c - x& ' for /=E CaXa ( 10 ) 
(a, v) = m 

where m = min{ (a, v) | a € A }. For the connections with Grobner bases, see [20] . For 
/ supported by A, denote the support of in v (/) by in v (A). A system / = (fi, f 2 , ■ ■ ■ , f n ) is 
supported on (A±, A2, ■ ■ ■ , A n ). 

We look for v so that in v (/)(x) = has solutions in (C*) n . A nonzero vector v is a pretropism 
for the system /(x) = if #dn v (/fc) > 2 for all k ranging from 1 to n. Every tropism is 
a pretropism, but not every pretropism is a tropism, as pretropisms depend only on supports 
A = (A 1 ,A 2 ,...,A n ) of/. 

The Cayley trick formulates a resultant as a discriminant as in [10\ Proposition 1.7, page 274]. 
Here we follow the geometric description of [19], see also §9.2]. Via the Cayley embedding we 
reduce A to one set: 

E A = {A x x {0}) U {A 2 x {ei}) U • • • U (A n x {e n _i}) (11) 

where is the k-th. (n — l)-dimensional unit vector. We call the convex hull of the Cayley 
embedding the Cayley polytope coiiv(Ea) 

We claim that enumerating all facet normals to coiiv(Ea) yields all tropisms. Furthermore: 
tropisms for curves are normals to facets spanned by at least two points of each support. In case 
of surfaces of solutions, we will find cones spanned by tropisms for curves. 

Algorithms to compute tropical varieties are described in [5] and implemented in Gfan [11] . 

Following from the second theorem of Bernshtein [3], the Newton polytopes may be in general 
position and no normals to at least one edge of every Newton polytope exists. In that case, there 
does not exist a positive dimensional solution set either. 

Running cddlib [9] to compute the H-repesentation of the Cayley polytope of the cyclic 8- 
roots problem yields 94 pretropisms. With symmetry we have 11 generators, displayed inTable[TJ 

For the cyclic 9-roots problem, the computation of the facets of the Cayley polytope yield 276 
pretropisms, with 17 generators: (—2, 1, 1, —2, 1, 1, —2, 1, 1), (—1, —1, 2, —1, —1, 2, —1, — 1, 
2), (-1, 0, 0, 0, 0, 1, -1, 1, 0), (-1, 0, 0, 0, 0, 1, 0, -1, 1), (-1, 0, 0, 0, 1, -1, 0, 1, 0), (-1, 0, 
0, 0, 1, -1, 1, 0, 0), (-1, 0, 0, 0, 1, 0, -1, 0, 1), (-1, 0, 0, 0, 1, 0, -1, 1, 0), (-1, 0, 0, 0, 1, 0, 0, 
-1, 1), (-1, 0, 0, 1, -1, 0, 1, -1, 1), (-1, 0, 0, 1, -1, 0, 1, 0, 0), (-1, 0, 0, 1, -1, 1, -1, 0, 1), 
(-1, 0, 0, 1, -1, 1, -1, 1, 0), (-1, 0, 0, 1, -1, 1, 0, -1, 1), (-1, 0, 0, 1, 0, -1, 1, -1, 1), (-1, 0, 
0, 1, 0, 0, -1, 0, 1), and (-1, 0, 1, -1, 1, -1, 0, 1, 0). 

The computations for n = 8 and n = 9 finished in less than a second on one core of a 3.07Ghz 
Linux computer with 4Gb RAM. For the cyclic 12-roots problem, cddlib needed about a week 
to compute the 907,923 facets normals of the Cayley polytope. Note that only a small fragment 
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pretropism v 



^solutions of in v (Cs)(z) 



-3,1,1,1,-3,1,1,1) 

-1,-1,-1,3,-1,-1,-1,3) 

-1,-1,1,1,-1,-1,1,1) 



94 
115 
112 
30 
23 
32 
40 
16 
39 
23 
509 



-1,0,0,0,1,-1,1,0) 

-1,0,0,0,1,0,-1,1) 

-1,0,0,1,-1,1,0,0) 

-1,0,0,1,0,-1,1,0) 

-1,0,0,1,0,0,-1,1) 

-1,0,1,-1,1,-1,1,0) 

-1,0,1,0,-1,1,-1,1) 



-1,1,-1,1,-1,1,-1,1) 



Table 1: Eleven pretropism generators of the cyclic 8-root problem, along with the number of 
solutions of the corresponding initial form systems. 



of the total number of facets yield pretropisms, so specific polyhedral algorithms for pretropisms 
will perform much better. 

5 The Second Term of a Puiseux Series 

In exceptional cases like the cyclic 4-roots problem where the first term of the series gives an exact 
solution or when we encounter solution lines like with the first four equations of cyclic 5-roots, we 
do not have to look for a second term of a series. In general, a pretropism v becomes a tropism if 
there is a Puiseux series with leading powers equal to v. The leading coefficients of the series is 
a solution in C* of the initial form system in v (/)(x) = 0. We solve the initial form systems with 
PHCpack [21]. For the computations of the series we use Sage [18] . 

In our approach, the calculation of the second term in the Puiseux series is critical to decide 
whether a solution of an initial form system corresponds to an isolated solution at infinity of the 
original system, or whether it constitutes the beginning of a space curve. For sparse systems, we 
may not assume that the second term of the series is linear in t. Trying consecutive powers of t 
will be wasteful for high degree second terms of particular systems. In this section we explain 
our algorithm to compute the second term in the Puiseux series. 

A unimodular coordinate transformation x = z AI with M having as first row the vector v 
turns the initial form system in v (/)(x) = into in ei (/)(z) = where ei = (1,0, ... ,0) equals 
the first standard basis vector. When v has negative components, solutions of in v (/)(x) = that 
are at infinity (in the ordinary sense of having components equal to oo) are turned into solutions 



Proposition 5.1. If the initial root does not satisfy the entire transformed polynomial system, 
then there must be at least one nonzero constant exponent ai forming monomial Cit ai . 

Proof. Let z = {z\, Z2, ■ ■ ■ , z n ) and z = (z%,Z3, . . . ,z n ) denote variables after the unimodular 



in (CT of in ei (/)(z) =0. 
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transformation. Let {z\ = t, Z2 = ri-, ■ ■ ■ , z n = r n ) be a regular solution at infinity and t the free 
variable. 

The ith equation of the original system after the unimodular coordinate transformation has 
the form 

/i = ^(Pi(z) + 0(zi)Qi(z)), i = l,2,...,n, (12) 

where the polynomial Pi(z) consists of all monomials which form the initial form component of 
fi and Qi(z) is a polynomial consisting of all remaining monomials of /j. 

For the second term in the series expansions, we denote 

Zl Z l w ■_ (13) 

Zi — T{ -j- k{t , 2 — 2, . . . , Tl. 

We first show that polynomial z™ l Pj(z) cannot contain a monomial of the form Cit ai on 
substitution of f j 1 3 p . The polynomial z^Piiz) is the initial form of fi, hence solution at infinity 
(zi = t,Z2 = T2,zz = r3, . . . , z n = r n ) satisfies z^P^z) entirely. Substituting (JT3J) into z^P^z) 
eliminates all constants in t mi Pi{z). Hence, the polynomial Pi(t) = Ri(t w ) and, therefore, t mi Pi{t) 
= Ri(t w+mi ). 

We next show that polynomial Qi(z) contains a monomial Cit ai . The polynomial Qi(z) can 
be rewritten as 

^Qi(z) = ^Tio(z) + z^Taiz) + • • • + z\^ n T m (z). (14) 

The polynomial Qi(z) = z?*Qj(z) consists of monomials which are not part of the initial form 
of fi. Hence, on substitution of solution at infinity (fT3)) . z\ l Q; L {z) = t ai Qi(t) does not vanish 
entirely and it must contain monomials which are constants. Since Qi(t) contains monomials 
which are constants, t ai Qi(t) must contain a monomial of the form Cijt ai . □ 

If the initial root does not satisfy the original system, then we may have a second term in the 
Puiseux series. Assume the following general form of the series: 

\ Zl (o) , . w ■ ( 15 ) 
[ Zi = c\' + yit w , t = 2,...,n, 

where c[°^ G C* are the coordinates of the initial root, is the unknown coefficient of the second 
term t w , w > 0. Note that only for some yi nonzero values may exist. We are looking for the 
smallest w for which the linear system in the j/j's admits a solution with at least one nonzero 
coordinate. Substituting (fT5|) gives equations of the form 



c?V*(l + O(t)) + t w+bi ^7y%(l + 0(t)) = 0, i = 1,2, . . . ,n, (16) 

3=1 

for constant exponents Oj, 6, and constant coefficients c| ^ and 7^. 

In the equations of (|16p we truncate the O(t) terms and retain those equations with the 
smallest value of the exponents because with the second term of the series solution we want 
to elimate the lowest powers of t when we plug in the first two terms of the series in the system. 



7 



This gives a condition on the value w of the unknown exponent of t in the second term. If there 
is no value for w so that we can match with w + bi the minimal value of ctj for all equations where 
the same minimal value of en occurs, then there does not exist a second term and hence no space 
curve. Otherwise, with the matching value for w we obtain a linear system in the unknown y 
variables. If a solution to this linear system exists with at least one nonzero coordinate, then we 
have found a second term, otherwise, there is no space curve. 



6 Series of Cyclic 8 and 12-Roots 



6.1 cyclic 8-roots 

We illustrate our approach on the cyclic 8-roots problem, denoted by Cs(x) = 
pretropism v = (1, —1, 0, 1, 0, 0, —1, 0). The corresponding initial form system is 



and take as 



in v (C 8 )(x) 



(17) 



X\ + X 6 = 

X\X2 + X 5 X 6 + X 6 X 7 = 

X4X5XQ + X5XQX7 = 

XqXiXqXj + X4X5XQX7 = 
XqXiX2XqXj + X0X1X5XQX7 = 

X0X1X2X5X6X7 + 2:03:1X4X5X6X7 + X1X2X3X4X5X6 = 
X0X1X2X4X5X6X7 + X1X2X3X4X5X6X7 = 
X0X1X2X3X4X5X6X7 — 1 = 

A unimodular matrix M £ ^ nxn j s an invertible integer matrix. We use unimodular matrices for 
coordinate transformations. Our unimodular matrices have the tropism as their first row: 



M 



-10 10 0-10" 



1 












10 

10 

10 

1 













1 





(18) 



Then the corresponding unimodular coordinate transformation, denoted as x = z M , is defined as 

X = Z , Xi = z 1 /z , x 2 = z 2 , X 3 = Z Q Z 3 , X 4 = Z4, x 5 = z 5 , x 6 = zq/z , x 7 = z 7 (19) 
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Applying x = z M to the initial form system ()1T|) gives 



in v (C 8 )(x = z M )= { 



(20) 



Z\ + Zq = 

ZiZ 2 + z 5 z 6 + Z 6 Z 7 = 
Z4,Z 5 Zq + 2 5 Z 6 2 7 = 
Z 4 Z 5 ZqZ 7 + Z\ZqZ 7 = 
ZiZ 2 Z 6 Z 7 + Z\Z 5 ZqZ 7 = 

Z1Z2Z3Z4Z5Z6 + ZlZ 2 Z5Z6Z 7 + ZiZ4Z 5 ZqZ 7 = 
Z1Z2Z3Z4Z5ZQZ7 + -21^2^4 Z^Z-j = 

. z\ZiZzZ 4 z^z%z 7 -1 = 

By construction of M, observe that all polynomials have the same power of zq, so zq can be 
factored out. Removing zq from the initial form system, we find a solution 



Z =t, Z\ = -I, z 2 



1 



I 

2' 



Z3 



-1, Z 4 = 1 + J, Z 5 



1 J 

2 + 2' 



z 6 = J, z 7 



-1 



(21) 



where / = \f— 1. This solution is a regular solution, i.e.: it has multiplicity one. We set z$ = t, 
where t is the variable for the Puiseux series. In the computation of the second term, we assume 
the Puiseux series of the form 



zq 


= t 


Z\ 


= -I + cit 


Z2 




Z3 


= -1 + c 3 t 


z 4 


= 1 + I + c 4 t 


Z5 


= \ + 1 + c 5 t 


z% 


= I + c 6 t 




= (-l-I) + c 7 t 



(22) 



We first transform the cyclic 8-roots system Cg(x) = using the coordinate transformation 
given by (|19p and then substitute the assumed series form into this new system. Since the next 
term in the series is of form Cji , we collect all the coefficients of t 1 and solve the linear system 
of equations. 

The second term in the Puiseux series expansion for the cyclic 8-root system, in the particular 
direction v, has the form 



zo 


= t 


Z\ 


= -/+(-l- 


Z2 


2 2 T 2 b 


Z3 


= -1 


ZA 


=1+1 -t 


Z5 


2^2 2 L 


Z6 


= I+(1 + I)t 


i Z7 


= (-l-I) + t 



I)t 



(23) 
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Because of the regularity of the solution of the initial form system and the second term of the 
Puiseux series, we have a symbolic-numeric representation of a quadratic solution curve. 

If we place the same pretropism on another row in the unimodular matrix, then we can develop 
the same curve starting at a different coordinate plane. This move is useful if the solution curve 
would not be in general position with respect to the first coordinate plane. 

For symmetric polynomial systems, we apply the permutations to the pretropism, the initial 
form systems, and its solutions to find Puiseux series for different solution curves, related to the 
generating pretropism by symmetry. 

In certain instances, the first term in the Puiseux series satisfies the entire system. We 
illustrate this situation by taking the pretropisms v = (1,-1,1,— 1,1,— 1,1,— 1) of the cyclic 
8-roots system and compute its initial form: 

x\ + x 3 + x 5 + x 7 = 

XqX\ + X X 7 + X\X2 + X 2 X 3 + X3X4 + X4X5 + X 5 X 6 + X 6 X 7 = 
X XiX 7 + X1X2X3 + X3X4X5 + X 5 X 6 X 7 = 
X XiX 2 X 3 + X XiX 2 X 7 + X XiX 6 X 7 + XoX 5 X 6 X 7 
+ X1X2X3X4 + X2X3X4X5 + X3X4X5X6 + X4X5X6X 7 = 



in v (C 8 )(x) = < 



1 •^i-^2^3^4 ' ^z-^3-^4^5 ' ■^i^l^b^o ~ ^4-^5 ^d-^Y — u 
XQX1X2X3XJ + XoXiX5X6X 7 + X1X2X3X4X5 + X3X4X5X6X 7 = 
X0X1X2X3X4X5 + X XiX2X 3 X4X 7 + X XiX 2 X3X 6 X 7 + XoXiX 2 X 5 X 6 X7 
+ XoXiX4X5X6X 7 + XoX3X4X5X6X 7 + X1X2X3X4X5X6 + X2X3X4X5XQX7 = 
X0X1X2X3X4X5X7 + XoXiX2X3X5X6X 7 + XqX\XzX^X^X%X 7 + X1X2X3X4X5X6X7 

X0X1X2X3X ^4X5X6X7 — 1 = 



= 



X0X1X2X3X4X5X6X7 

The coordinate transformation for the initial form system 
matrix 



(24) 

is given by the unimodular 



M 



-1 
1 









The coordinate transformation x - 

z 3 /z , x 4 = z Zi, x 5 = 2:5/^0, x 6 = z z e 



1 



1 







M 



-1 




1 








1 





1 







-1 






1 






1 







1 





-r 








1 



yields x = z , x x = Zi/z , x 2 = z z 2 , x 3 = 
X7 = Zj/zq. Applying the coordinate transforma- 
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tion to the initial form system (|24p gives 



in v (C 8 )(x = z 



M x 



Z\ + z 3 + z 5 + z 7 = 

Z1Z2 + Z2^3 + Z 3 Z 4 + ZAZ 5 + Z 5 Z 6 + Z 6 Z 7 + 2i + Z 7 = 
21^2^3 + Z 3 Z4,Z 5 + Z 5 Z 6 ^7 + ZlZ7 = 
ZIZ2Z3Z4 + Z 2 Z 3 ZiZ 5 + 2 3 Z 4 Z 5 Z 6 + Z 4 ^5^6^7 + ZiZ 2 Z 3 
+ ZiZ 2 Z 7 + ZiZqZ 7 + Z 5 ZqZ 7 = 

Z\Z 2 Z 3 Z 4 Z^ + Z 3 Z 4 Z 5 Z 6 Z 7 + ZlZ 2 ^3^7 + ZiZ 5 Z 6 Z 7 = 
^12 2 ^3^4^5^6 + Z 2 Z 3 Z i Z 5 Z & Z 7 + ^1^2^324^5 + Z X Z 2 Z 3 Z^Z 7 
+ ZiZ 2 Z3Z6^7 + 21 Z 2 Z 5 Z$Z 7 + ZiZ 4 Z 5 Z e Z 7 + 23Z4Z 5 Z 6 Z 7 = 
ZlZ 2 Z 3 Z 4 Z 5 Z 6 Z 7 + Z1Z2Z3Z4Z5Z7 + Z 1 Z 2 Z 3 Z 5 ZqZ 7 + Z 1 Z 3 Z 4 Z 5 Z 6 Z 7 

ziz 2 z 3 z 4 z 5 z 6 z 7 -1 = 



(25) 



The initial form system (125p has 72 solutions. Among the 72 solutions, a solution of the form 
z =t, zi = -1, z 2 = I, z 3 = -I, Zi = -1, z 5 = 1, z 6 = -I, z 7 = I, (26) 
here expressed in the original coordinates, 

Xo = t, X\ = —1/t, x 2 = It, x 3 = —I/t, X4 = — t, X5 = 1/t, Xq = —It, x 7 = I jt (27) 

satisfies the cyclic 8-roots entirely. Applying the cyclic permutation of this solution set we 
can obtain the remaining 7 solution sets, which also satisfy the cyclic 8-roots system. 

In |23j . a formula for the degree of the curve was derived, based on the coordinates of the 
tropism and the number of initial roots for the same tropism, and is given in the definition below. 

Definition 6.1 (Branch Degree). Let v = {v\,v 2 , • • • ,v n ) be a tropism and let R be the set of 
initial roots of the initial form system in v (/)(z). Then the degree of the branch is 



max Vi 



rt 

rn.va.Vi 
i=l 



(28) 



We apply Definition 16.11 to the Puiseux series developments for the cyclic 8-roots problem and 
obtain 144 as the known degree of the space curve of the one dimensional solution set. 

(1,-1,1,-1,1,-1,1,-1) 8x2 = 16 

(1,-1,0,1,0,0,-1,0) (1,0,0,-1,0,1,-1,0) 8x2 + 8x2 = 32 
(1,0,-1,0,0,1,0,-1) (1,0,-1,1,0,-1,0,0) 8x2 + 8x2 = 32 
(1,0,-1,1,0,-1,0,0) -> (1,0,-1,0,0,1,0,-1) 8x2 + 8x2 = 32 
(1,0,0,-1,0,1,-1,0) -> (1,-1,0,1,0,0,-1,0) 8x2 + 8x2 = 32 

TOTAL = 144 



Table 2: Tropisms, cyclic permutations, and degrees for the cyclic 8 solution curve. 



Using the same polyhedral method we can find all the isolated solutions of the cyclic 8-roots 
system. 
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In the next few paragraphs, we address, informally, some of the issues related to the time 
complexity of the procedure for computation of the second term. 

In the direction of tropism (1, —1,0, 1,0,0, —1,0), there exists a second term in the Puiseux 
series. The procedure solves the initial form system, which yields 40 solutions, and checks whether 
the first term satisfies the cyclic 8-roots system. Then it proceeds to construct and compute the 
second term in the Puiseux series. The total time required is 35.5 seconds, which includes 28 
miliseconds that PHCpack needed to solve the initial form system. 

For the tropism (1, —1, 1, —1, 1, —1, 1, —1) there is no second term in the Puiseux series. In 
this direction, the first term solves the entire cyclic 8-roots system. Hence, the procedure for 
construction and computation of the second term does not run. It takes PHCpack 12 seconds to 
solve the initial form system, whose solution set consists of 509 solutions. Determining that there 
is no second term for the 509 solutions, takes 199 seconds. 

Given their numbers of solutions, the ratio for time comparison is given by ^ = 12.725. 
However, given that for tropisms (1, —1, 0, 1, 0, 0, —1, 0) the procedure for construction and com- 
putation of the second term does run, unlike for tropism (1, —1, 1, —1, 1, —1, 1, —1), the ratio for 
time comparison is not precise enough. A more accurate ratio for comparison is ~ 5.686. 

6.2 cyclic 12-roots 

The generating solutions to the quadratic space curve solutions of the cyclic 12-roots problem are 
in Table [3l As the result in the Table [3] is given in the transformed coordinates, we return the 
solutions to the original coordinates. For any solution generator (n, r 2 , • • • , rn) in Table O 

z = t, zx = n, z 2 = r 2 , z 3 = r 3 , z 4 = r 4 , z 5 = r 5 , 

(29) 

zq = r 6 , z-j = r 7 , z 8 = r 8 , z 9 = r 9 , z 10 = r w , z n = r n 
and turning to the original coordinates we obtain 

x =t, xi = n/t, x 2 = r 2 t, x 3 = r 3 /t, x 4 = r 4 i, x 5 = r 5 /t 

x 6 = r 6 t, x 7 = r 7 /t, x 8 = r 8 t, x 9 = r 9 /t, x 10 = r 10 t, x u = m/t 

Applying Definition 16. 1\ we see that all space curves are quadrics. Compared to [16] , we arrive 
at this result without the application of any factorization methods. 



7 Concluding Remarks 

Inspired by an effective proof of the fundamental theorem of tropical algebraic geometry, we 
outlined in this paper a new polynomial method to compute Puiseux series expansions for solution 
curves of polynomial systems. The main advantage of the new approach is the capability to exploit 
permutation symmetry. 

For our current preliminary implementation we relied on cddlib for the pretropisms, the 
blackbox solver of PHCpack for solving the initial form systems, and Sage for the manipulations 
of the Puiseux series. Because of the many spurious normals caused by the Cayley embedding, 
the pretropism calculation can still be improved significantly. Often it may happen that the 



12 



initial form systems are sparser than the original system and are thus easier to solve. Another 
improvement to the new approach would be to replace a blackbox solver by a recursive application 
of the new approach to the initial form systems. 
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Table 3: Generators of the roots of the initial form system in v (Ci2)(x) = with the tropism v 
(+1, — 1, +1, — 1, +1, — 1, +1, — 1, +1, — 1, +1, — 1) in the transformed z coordinates. Every solution defines a solution curve 
the cyclic 12-roots system. 



